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One of th.e vital operating attitudes of an airplane 
is the time of take-off, and the question arises as to how 
the dimensions of the landing gear affect the stability of 
an airplane while rolling in horizontal position. 

The prohlems of the stability theory, presumedly known, 
refer to the free, level and three-dimensional motion with 
three and six degrees of freedom. Ordinarily, they are too 
complicated to afford immediate practical results for the 
designer; for example, a dynamical stability investigation 
rarely exerts any essential effect on the design of an air- 
plane . 

The stability analysis of an airplane while rolling is 
much more simplified to the extent that it can be obtained 
for niTmerical data which can be put to practical use in the 
design of landing gear dimensions. Every landing gear type 
attains to a critical ground friction coefficient that de- 
cides the beginning of instability, i.e., nosing over. 
This study has, in addition, a certain interest for the use 
of wheel brakes. 

Several premises are called for. The airplane rolls 
over flat, level ground and at steady speed u. The sta- 
bility of this motion is to be examined by adducing dis- 
turbances to this originally uniform motion and follow its 
course in time. The speed maximum u is assumed such that 
the gross weight is partly carried by wing lift, and partly 
supported by the ground. A second premise is that only the 
wheels touch the ground and the longitudinal airplane axis 
is parallel with the ground surface. 

While rolling the airplane executes a mechanical mo- 
tion with two degrees of freedom, which can be forced by a 
certain elevator and throttle setting. It. is clear that a 
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comprehensive analysis of the stability wo-uld have to in- 
clude the actually accelerated and- decelerated motion at 
take-off and landing; hxit mathemat icail difficulties pre- 
vented this general invest igation from being completed. 
The shock absorption was disregarded fOT similar reasons. 

Another important problem is the stability at landing 
when landing gear and tail skid already touch the ground. 
This study, Significant for brake design, can be reduced 
to the relative balance of the forces, inertia included, 
acting on the airplane; but this is beyond the scope of 
the present repor'^. 

Boiling of the airplane with constant speed u . is 
expressed by the fpllqwing equations of equilibrium; 



.Herein- 3 = thrust, (J = gross weight, N = wheel 
pressure, R = friction, p = air density. Cat c'^ and 
Cjn the well-known air coefficients, P = wing area, and 
I = distance of center of pressure 0 from the axis of 
the tail surfaces; a' and b' denote rear and height po- 
sition of point of contact B of the wheel with reference 
to a system of coordinates, the origin of which is coinci- 
dent with the e.g. of the airplane. 

". The f r ict ion is 



f being the friction coefficient. 

As to the prefixes used in equations 1-3, be it 
noted that the x direction in the direction of motion, 
the y direction is positive downward. Clockwise moments 
are positive. Cjj^ positive signifies a nose-heavy moment 
and carries - (minus) sign. 




(3) 



(2) 



(1) 



E = f N, 



The disturbed motion of the rolling airplane is rep- 
resented as plane motion by statement of the speed of an 
airplane point and by the angular velocity around this • 
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point, i.e., abo-al; the axis which passes through this 

point and is at right angle to the plane of motion.' Ad- 
mittedly, the speed changes as, the, .reference point changes, 
whereas the angular velocity remains constant.' 

Disregarding the rolling of the wheels, the velocity 
component of contact point B is horizontal. The motion 
of the airplane can he expressed "by the horizontal, velocity 
of point B and hy the angular velocity ahout this point. 



As the airplane changes its originally horizontal po- 
sition through an angle -B- = 0 in smooth rolling), the 



angular 

the e.g. 
and Vg, 



d^ 
= dt' 



velocity is 

, has the horizontal 
interdependent with 



( t time). 

and 
u 



vertical 
and ^ . 



u. 



= u - b' 



dt 



The. velocity of 
components u- 



(4) 



dj^ 
dt 



= positive for positive rotation of airplane, that is, 
fore part of airplane rises.) 

The horizontal rolling motion is disturbed by the 

change in horizontal speed ^u of point B and by the 
angle of slope The angle of attack changes as Aa, 

wher e 



A- = * + 5' I? 



(5) 



The change in angle of attack, induced by the turning 

motion of the airplane, is approximately put at — ' 

u dt 

The equations of motion of the airplane can be ex- 
pressed in the previously introduced system of coordinates 
as follows: 



du 



m 



m 



s 



dt 
dv 

s 

dt 



2 

J — _ = M-n«V 



dt 



dt 



= s - E - AE - w + A a' |i 

u dt 

= G - N - AN - A -. W a' |# " 

. , . .. u dt 

•+••: (N+AlT) (a'+b« ^) - (E+AE) (b«-a«^) 



(6) 
(8) 
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Herein 
to = ma s s 



J ■■ 



motoent i&f- infertia 



of airplane 



,A =. lift 
117 = dras 



ch&ngeiS in, wheel pressure and friction during 
di-s-turhed motion 



M 
V 
n« 



moment qf. a,ir loads referred to' e.g. 
resultant airplane speed' 
damping coefficient 



In connection with these equations, lift and drag form 
a2igi;e ^ - <9 , and cp with the horizontal, wherehy 

qp = arc tan — — 
u 

s 

in addition, we assumed cos cp = 1 and sin 0) = — '~ . 

^ ^ ■ u dt 

(The relative speed is sho^rn in PigurS 2.) Due alloisrance 

was made in the moment equation, that airplane rotation d- 

modifies lengths a" and h' to a' + h'^ and "b' - a>^. 



Furthermore, V, A, W, and M are expressed as 
= u^- 3 u h' If + (a'2+ ^•=)(|f) 



— 2 2 , 2 

V = Ug + Vg 



dt 



A = b; p V (a + Aa) 



2 - ■ -a 

3 



(9) 



W = I ! V c^ (a + Aa) 

K = I P V^Cjn (a + Aa) I 

where a = an^le of attack in equilihrium attitude. 

Coefficients (a + Aa) etc., are written Cg^ + 
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-r-=- Aa, etc., and -3 — etc., later designated "by Co etc. 
da da » a 

Kote, that A ti, , M are; small qtiantitiea whose 

dt 

powers and products can he disregarded. With a thought 
to (1), (2), and (3) equations (6), (7), and (8) now he- 
come . ■ 



2 

+ Au p P u c^ + f AN - h» m 



at 
dt 



2-2 ^°w3 



+ ^ I P u^ 6^ = 0 



Au pP u Ca + AN + a'm 



d=i^ 



dt 



dt s 

pP u[— 5— + — — - t c^] 



(10) 



+ ,^ §■ F u'' c= = 0 



Au pP u I +Alir (f b» - a») + J 



dt^ 



a ' c 



2^ - ^'°m + 



dt ^ ^ t 

+ [ I P Cjj^ - N (a'f + h»)] = 0 



(11) 



(12) 



Here it vtas already taken into account that AE = f AN. 
There are two degrees of freedom for this guided airplane 
motion, whereas it requires three equations hecause of the 
appearance of the unknown force of reaction AN. 



As known, the solution for Au, AN and S- is written 
in form K e'^'' , wherein t = time and K and K are con- 
stants. The definition of \ yields. the algehraic equa- 
tion 

inX+pPu c^, f h»m+\pPu [- ^^+^^ -'^'0^ ]+ | Pti? c. 

a • a » c. 



pPufe^; • .1 +\2a'mt\ pP u [~^+-2-^ -^'°a 3 + 2 



pPu I Cjjj (fh«-a«) \2 j+x pi.^a[^™-h'cm+n«3 + 



+ P -0-2 

. 2- 



Pu^ I - N (a' f + h) 



= 0 



(13) 
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SubstitTiting r\ for X , where 



•n = —5— \ 

equation (13) assumes, after several tfearrangement s , the 
form of 



•n + 2c^ f 
2c„ 1 



^ ® " °a^ + C 

+ Ti a + c^) + fia ° 



2c 



m 



f h - a Tf D + TiCa&jn + 2n] 



N 



m 



pl-a 2 



L_[af + -b] 



= 0 



(14) 



Here we put 
I 



-a' T, - « - 5l' p - -! 



2m_ 
pF I 



m I' 



Quantities c„ and -—5- can he defined from (2) 

p Pu"* 

and (3). With the introduction of the nondimensional 
quantity = . 



Cjn = (a - f t) (gg - c^) 



pPu" 



So - 



(15) 



is bhtained. 



Whan these quantities are Inserted in determinant 
(14), the latter becomes 

cp^ + Ti^ + 9^ T) + q)^ = 0 (16) 

Coefficients <5P * <P > 9 ^ cp can he written as: 

0 1 33 

= a^ - ah f + D 



a^[f 2 (go - da) + + Sc^] - ah [ f 2 (go - c^) 
■ + f (Ca + 3c^)] + a + 2n + D(2c^ - 2caf) 
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= a^Cf (2c^ - 2c^c^ + ScaSo - 

+ 2c^3 + 2S^c^ + 2gp ( - c^)] 

- a b [2 f=(go Ca) (c^ + c^) 

+ f (ac^^ + 2c„c ) + 2 g„ (c„ - c„)] 
a a w 0 a W-" 

+ a CSCjjjC^ - 2f c„cj + 4 c„ n - 4 f c^n 



'm^a 



+ C(aCc^ - fg^ + fcj + b Cc^ - g^ - c^f] + cj 
2C [a(-f^c^ (g^ - Ca) + fCg^ - Ca) (c^ - c^) 
+ ^a°w - ^Sq) - ^(^^°a^So - c^) 
+ fCc^c^ - So °a <So 
+ ^mCc^f + V)] 



(17) 



The condition for the stahility is, that the inequa- 



t i ons 



%■ > 0, cp^ > 0, cp^ > 0, cp_^ > 0 
are complied with. 



(18) 



The significance of these stability criteria can also 
he so expressed that the aperiodic proportion' of the- dis- 
turbed motion changes from damping to intensification, when 
the condition 0 is complied with. This condition is 

called static stability; the root in (l5) is real, the 
transition frbm damping to intensification is noted as 
T\ = 0 . 

The condition ^j_^>2 - 'Pq'^'s ~ ^ equivalent in so 

far as T\ becomes purely imaginar jf- ; i . e . , ■ the periodic 
portion of the disturbed motion is neither damped nor in- 
tensified. 

The primary . aim of thes e . stgbll ity . invest igati ons is 
to ascertain how the ■ stability of the rolling airplane is 
influenced by the dimensions of the landing gear. The aero- 
dynamic factors, the constants C and D depend on the . 
general character Of the airplane. There remain then the 
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quantities a and b which denote the arrangement of 
the landing gear, and lastly, the coefficient of ground 
friction f. 

The stability conditions can be interpreted as func- 
tions of quantities , a, b and f. These functions can be 
utilized in two directions: one, to examine individual 
airplanes with definite a and b values and assuming 
the friction coefficient as variable. Then the latter can 
be numerically .defined from the stability equations. A 
condition of the ground conformably to this coefficient, 
produces instability when rolling. This may be called 
"critical friction coefficient," and with this factor, 
different airplanes, i.e., types of landing gears .can be 
compared. 

A second application of the stability investigation 
is a more general study of the relationship between a, b, 
and f in the manner that the effect of size of the land- 
ing gear on the stability of rolling motion is to be prin- 
cipally explained. 

Apart from the cited structural quantities a, b, and 
f, two parameters, c^^ and are equally important. The 

first is the criterion of the static stability of the air- 
plane in flight, the second depends on the speed of the 
rolling airplane. 

The equation of static stability of rolling motion of 
the. airplane, according to (17) and (18) is: 

^ = a C-f^ Ca(go - ca) + f (go - (c^ - c^) 

+- °a% - Vo ^" ^ [^^^°a(So - °a> 
+ f(c^c^ - C^go + Ca(go - Ca^> + ^w^So - °a^l 

+ 6j„(c^f + c^) ^0 (19) 

An analysis of this (19) reveals the startling effect 
of a, b, f and c-^. When a > 0, that is, when the e.g., 
is back of the landing gear, the static stability is in- 
creased, when b becomes large, i.e., the landing gear 
too high, the static stability decreases. Static stability 
is also increased when the airplane is stable in the usual 
sense of the work, i.e., when c >0. 
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It sh.oTild te remembered that a may "be increased on- 
ly- to,..a certain degree in the positive sense, or else it 
might happen that the tail of the airplane cannot ..be raised 
with the elevator at the right time, when taking off. 
Q,xiantity b is bo-anded by the propeller. 

The Tinstabil izing; action of the friction becomes read- 
ily apparent. . An increase in angle of attack increases the 
lift and decrease.s the .wheel pressure and the friction. 
The decrease- in. friction denotes a tail-heavy moment in the 
same, direction as the direction of rotation of the airplane. 

. These problems were studied by means of concrete exam- 
pl,es. Part of the data is given in Figures 3 to 6, 

The airplane' characteristics were: 

G = 1,500 kg, J = 300 kgm/s^, P = 22m^, I = 4.9 m 

C = §-S. = 22.6, ■ D ~. -^^= 0.08 
pFl m 

c^ = 0.43, Qa =.3.9, c^ = 0.06, c^ = 0.200, 

•c = 6.128, . n = 1 .09 . 
m 

Figure 3 shows function 'Pg = 0 in the a, b plane; 
f is considered as parameter and g^ = 0.63. The influ- 
ence of the friction is manifest. To examine a concrete 
landing gear type, use the critical coefficient of fric- 
tion of Figure 3. At low friction the landing geal* could 
be fitted aft of the e.g., but as the friction increases 
such an arrangement becomes impractical. The influence of 
the forward position of the lending gear can be seen in 
Figure 4, where the increase in critical coefficient of 
friction by increasing' a is distinctly noticeable. To 
examine the influence of the rate of rolling the stability 
equation ^3. = 0 is pictured for gjj = 0.50, i.e., at a 

higher speed. As is seen in Figure 5, the critical fric- 
tion numbers are higher for g^ = 0.50 than for g =0.63 
(see fig. 3) with fixed a, b values. The static stabil- 
ity grows with the speed. 

In conclusion the influence of c^j^ of the sta"tic "sta- 
bility in flight is to be examined. For comparison consid- 
er a neutral airplane that differs from the one treated by 
the value Cm = 0. According to Figure 6 the straights for 
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the different friction values meet in point a = b = 0, 
that is, they have shifted parallel in comparison to the 
stahle airplane. The critic,ai friction factors for the 
fixed' a, h values have hecbme' smaller; the stability in 
rolling is lower.' 

By analyzing the periodic solutions of the motion 
equations, it is found that stability prevails when the 
periodic parts have a damping factor. The transition from 
damping to intensification occurs -when the root of the 
typical equation (16) becomes purely imaginary, the condi- 
tion o.f which is: . 



= 0 



(20) 



An approximation is to afford a survey relative to 
this condition and its importance from the des igners ' point 
of view. 



Function (Pg can be simplifiesd by a consideration of 
arrangement of quantities. The terms of the function mul- 
tiplied by 0 are large compared to the others; for in- 
stance, put t ing quant it ies a, b and f, which practically 
range between 0 and 0.5, at =■ 0.3, the terms multiplied 
by C yield «- 19..5, the others - 0.202; accruing there- 
fore, when neglected, in a relative error of the order of 



10" . A satisfactory approximate calculation is 
a[ca - f (go - Ca)] + l^Cca - &o " °a^^ + ^ 



C 



m 



(21) 



analysis of 



Further simplifications may be attained to from an 

®ij£3_ and ^oILs. . Quantities a, b and f 
c c 

practically range between 0- and 0.5; calculation shows 
— to be of the order of size of 2 t o 3 between the 



cited practical limits, and 



^P^3_ to be of the order of 
C 

0.01 to 0.06. In the practical case, cp^^ is alv/ays posi- 
tive. 



To arrive at a summary regarding the zero points 
the stability equation (dynamic stability) 
the ifiore siiiiTDle 



cp CD 



of 



1b. 



(22) 
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•will' have to "be examined, which is decisive so that. (20) 
may "become = 0 in the practical range. The proof .of .this, 
stat'enient follows from the consideration of the ar.range- 
me'nt of the quaiiitit ies . Of' course, the acctiracy of this- 
approximation mtist be checked when the order of size of 
the numerical val-ues is at appreciahle variance with the 
erstwhile quoted data. The transition from dynamic sta- 
bility to instability is thus characterized by the- equa- 
tion: • 

9^ = 0 

This equation is illustrated like the criterion of 
the static stability. The coefficient of friction, consid- 
ered as parameter, is put at 0<f <0.5. Each parameter 
value has a corresponding linear function of a and ;b. 
These straights, plotted in Figure 7, readily show the de- 
pendence of the critical coefficient of friction on the di- 
mensions of the landing gear. 

As the coefficient of friction increases a consistent- 
ly higher a value is necessary to assure stability; i.e., 
the landing gear must be ahead of the e.g. The influence 
of it will be noted, is similar to that established by 

the examination of the static stability. 

Some remarks about are to form the conclusion of 

this report. By 6 was merely meant the measure of the 
static stability, the differential quotient of the moment 
curve, and in the numerical example this value was taken 
from the usual stability calculations. The process must 
be looked upon as an approximation, the accuracy of which 
becomes questionable, because of the motion of the airplane 
in proximity of the ground. This factor can be accounted 
for by reflecting the airplane respectively, the Prandtl 
vortices of the wings and of the horizontal surfaces about 
the plane of the ground and then effect the moment and 
downwash corrections. 

These calculations do not obtain to a complete de- 
scription of the process on the horizontal tail group. Be- 
cause of the slipstream a twofold action on the horizontal 
tail group is noticeable. The velocity of flow is in- 
creased; the induction effe(3t of the wings on the control 
surfaces is influenced by an area of discontinuity formed 
by the jet boundary. An exact recognition of the phenom- 
ena is moreover, rendered difficult because of the exis- 
tence of a radial motion in the slipstream. 
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The quot ed ^ influences should- really he amenahle to 
estimation If Cj^ were to he correctly introduced in the 
cal culat ions i Apar<ijt- from these refinements, the applica- 
tion of the stahility theory affords a correct numerical 
study of the interference of the landing gear designs and 
of the state or conditions of the ground. 

Translation hy J . Vanier, 
National Advisory Committee 
for Aeronautics. 



IT.A.C.A. Technical Memorandum No. 682 





N.A.C.A. Technical MemorandGun No. 682 



Figs. 3,4,5 




93[a,b,f]=0 
go = 0.63 

Pig. 3 function <P3 = 0, coefficient of friction f as parameter. 
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a 



■b= 0.326 

<+53ra,0.326,f] = 0 



Fig, 4 Effect of forwaxd position of landing gear. 




cpg [a,b,fl = 0 
go=0.50; 



Cjjj= +0.128 



Fig. 5 Function go=0, coefficient of friction f as parameter 
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cp [a,-b,f]= 0 
go =0.50 



m 



lig. 6 Neutral airplane, 0^= 0. Tlie straigh.ts meet in the 
origin of the coordinates, that is, are displaced 
parallel compared to the stahle airplane. 




cpg [a,h,f] = 0 



5ig, 7 Critical coefficient of friction plotted against 
laQ^lng. gear dimensions. 
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